Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex

Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.